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ABSTRACT 


In this work we mathematically reconstruct and document the digital position 
controller implemented in the control computer of the 3-axis attitude motion 
simulator, since the information supplied with the executable code of this controller 
was insufficient to make substantial modifications to it. We also develop 
methodologies to introduce changes in the controller which do not require rewriting 
the software. Finally, recommendations are made on possible improvement to the 
control system performance. 
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INTRODUCTION 


The objective of this project is to investigate the nature and 
details of the implementation of the digital position controller of the 
three-axes attitude motion simulator in order to effect changes in its 
parameters without having to rewrite the software, since the source 
code is not available. 

To get a good understanding of this digital position controller, 
several experiments were carried out on the simulator via the control 
computer, with the assistance of NASA/Boeing personnel. The goal of 
these experiments was to determine whether or not the position con- 
troller implemented in the control computer was the discrete-time 
equivalent of the position controller implemented for the analog mode 
operation of the simulator [1], 
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DETERMINATION OF THE DIGITAL POSITION CONTROLLER 


Daring the analog mode operation of the simulator, the position 
controller of any of the three axes has the form 


C fc (s) 


A U (s) 
E(s) 


K(s/oj^+ 1 ) (s/u^+l) 
s (s/w^+l ) 


( 1 ) 


where E(s) is the Laplace transform of the error signal e(t) and U(s) is 
the Laplace transform of the control signal u(t). The error signal e(t) 
is defined as the difference between the commanded angular position 0 (t) 
and the actual system angular position 0(t). 

To obtain a discrete- time equivalent of this analog controller, the 
industry-accepted way of doing it is by applying the bilinear transforma- 
tion, since it preserves the stability properties of its continuous-time 
ccun t erpar t , i . e . , 
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where 


E(z) = Z-transform of the discrete-t ime error signal e(kT), 

U (z) = Z-transform of the discrete-time control signal u(kT), 

A K W 3 

d ■ r. - : » 


( 3 ) 


XIX- 2 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(4) 


A 

a o 


(!) 2 


+ Y + a) i t0 2 


A 

a l = 


(!) ♦ ! »3 

{(f) 2 - -.-J . 


(!) - ! 


a 2 = 


A (|) - | ^i^ 2 } + V 

' (D 2 


, 2 

+ T W 3 


4 

K A T “3 
b l =— 5 


(!) 

(!) 2 - 


, 2 

+ T 3 


u A V T / T W 3 

b 2 * - 7^2 


(!) + ! 


0), 


T = sampling interval (in seconds). 


(5) 


( 6 ) 


(7) 


( 8 ) 


From equation (2) we find that the output of the digital position 
controller at the sampling instant kT is given by 


u(kT) = - bj u(kT-T) - b 2 u(kT-2T) + K^eCkT) + K^etkT-T) 
+ K Q a 2 e (kT-2T) 


(9) 


According to the information available on the implementation of the 
digital position control algorithm [1], the parameters of the controller 
can be changed from the keyboard when INIT (initialization) is selected 
from the control computer main menu screen. When this entry is selected, 
a total of nine constants (labeled "coef f icient s") per axis appear on 
the screen; however, it is known that only seven of them correspond to 
the digital position control algorithm, one to the digital rate command 
scaling, and one is set to zero. 
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As shown in equation (9), there is a total of five coefficients, 
which strongly suggests that the digital position controller implemented 
in the control computer is not the one obtained by applying the bilinear 
transform to the analog position controller transfer function given by 
equation ( 1 ) . 

It should be pointed out at this stage, nonetheless, that there are 
other realizations of the digital position controller given by equation 
(2) and that more constants could be generated for scaling purposes to 
avoid the usual problems encountered when fixed point arithmetic is used. 
However, the numerical values of the coefficients of these other realiza- 
tions do not match those that appear on the INIT screen for any of the 
three axes. 

To figure out the actual structure of the controller, several 
experiments were carried out to determine, for example, if integration 
of the error signal was taking place to make sure that the controller 
integrator had been implemented. Experimental results indicated that 
integration of the input signal did indeed occur, therefore, to generate 
the exact seven coefficients per axis, it became apparent that either a 
pole or a pole and a zero had been added to the original analog position 
controller . 

Let us modify the analog position controller by adding a pole at 
s = -o)^ , i. e . , 


(s/u>^ + 1) (s/u^ + 1) 
G e (s) ^ s(s/u)^ + 1) (s/to^ + 1) 


( 10 ) 


Then, application of the bilinear transform to this controller transfer 
function yields the following discrete-time equivalent controller 
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( 13 ) 
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A direct time domain realization of this digital controller is 
given by 
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u(kT) = - d lU (kT-T) - d 2 u(kT-2T) - d 3 u(kT-3T) + K^c Q e(kT) 
+ K^eCkT-T) + K^c 2 e(kT-2T) + K^e (kT-3T) 


( 20 ) 


This direct realization of the digital position controller is by 
no means the best realization as it is susceptible to round-off and 
quantization errors. However, the numerical values of its coefficients 
were found to be the same as those that appear on the INIT screen for 
each of the three axes once was determined, thus leading to the 

conclusion that the controller given by equations (11) and (20) is the 
one that has been implemented in the control computer (PC). In fact, 
the extra pole (-a> 4 ) is located at -20, -25, and -22.8 for the roll, 

yaw, and sidereal axis, respectively. 

The following analog transfer functions were determined to have 
been the source of the discrete-time, equivalent position controllers 
implemented in the control computer. 


Roll: 


G’(s) 


0 0047 (s/0.056 4- 1) (s/0.66 + 1) 
' s(s/5.6 + 1) (s/20 + 1) 


( 21 ) 


Yaw: 


G c (s) 


0 n?s (s/0.095 + 1) (s/0.588 + 1 ) 
s(s/l . 43 + 1) (s/25 + 1) 


( 22 ) 


Sidereal : 


G’(s) 


0 no (s/0.068 + 1) (s/12.5 + 1) 
s(s/20 + 1) (s/22.8 + 1) 


(23) 


Comparing equations (21) to (23) with those found on page 3-13 
of Reference 1, we can conclude that the analog position controllers 
were modified by adding a pole before they were discretized to obtain 
the digital position controllers. 

Application of the bilinear transform (with T = 0.01 sec) to equa- 
tions (21) to (23) yields the following discrete-time equivalent position 
controllers. 

Roll : 

C< ( z ) - 6 . 32005xl0~ 2 - 6 . 274936xl0~ 2 z~ 1 - 6 . 320026xl0~ 2 z~ 2 + 6 . 274959xl0~ 2 z~ 3 
1 - 2.763707 z -1 + 2.537319 z -2 - 7.736116 x 10 _1 z~ 3 


( 24 ) 
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Yaw: 


G c (z) 


7.08467x10 2 - 7.036407x10 2 z 1 - 7.084631x10 2 z 2 + 7.036446x10 2 z 3 


1 - 2.763579 z 1 + 2.530314 z 2 - 7.667345 x 10 1 z 3 


(25) 


Sidereal 


G'(z) = 


6.979685x10 2 - 6 . 153802xl0~ 2 z 1 - 6.979127x10 2 z 2 + 6.154359x10 2 z 3 


,-1-3 


1 - 2.613514 z 1 + 2.26424 z 2 - 6.507263x10 z 


(26) 
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Then the digital position controller for any of the three axes is given by 


G c (z) 


-1 -2 

K, + K 2 z + K 3 z + K^z 
] - K 5 z _1 - K^z ” 2 - K 7 z _ 


(34) 


The reason for writing the digital position controller in the form 
of equation (34) is that when the INIT menu selection is made, the 
coefficients of the digital position controller appear on the control 
computer monitor screen as through K^. For example, the coefficients 

of the digital position controller for the roll axis are: 
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Kj = 6.32005 x 10"* 

K 2 = -6.274936 x 10 _2 , 

K 3 = -6.320026 x 10~ 2 
= 6.274959 x 10 2 , 

K,. = 2.763707 

K 6 = -2.537319 

K ? = 7.736116 x 10 -1 

The present implementation of the digital position controller is 
therefore given by 


u(kT) - K^u (kT-T) + K 6 u(kT-2T) + K ? u(kT-3T) + K^eCkT) + l^e (kT-T) 


+ K 3 e(kT-2T) + K 4 e(kT-3T) 


(35) 


Remark : The dc gains of the controllers given by equations (21) 

through (23) are smaller than those of the original analog position 
controllers (shown on page 3-13 of Reference 1) because of the fact 
that the position sensors used in the analog and digital modes are 
different. In the analog mode, a single turn potentiometer is used to 
sense the angular position. This potentiometer has a gain embedded in 
it, in fact, it is equal to 0.05 Volts/degree for the roll axis and 
0.1 Volts/degree for the yaw and sidereal axis. In the digital mode, 
a resol ver- indue tosyn pair, along with two resolver to digital conver- 
ters, are used to obtain the angular displacement. The resolver to 
digital converters have a third-order dynamic model with unity dc gain. 
Therefore, this change in the gain is reflected by the reduction of 
the gain of the modified analog equivalent position controllers of 
equations (21) through (23). 
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CONCLUSIONS AND RECOMMENDATIONS 


The structure of the three-axes attitude motion simulator digital 
position controller implemented in the control computer has been deter- 
mined through experimentation and educated guesses, since no documenta- 
tion about (or computer program source code implementation of) such a 
controller is available. 

Thus, given that only the executable code of the controller is 
available, its structure cannot be changed without rewriting the entire 
control software; however, it is possible to make some on-line changes 
which do not require new controller software. For instance, the original 
location of the analog equivalent position controller poles and zeros 
could be changed to new desired values -w j , -u)£, -oo^, and and then 

compute through K j off-line using equations (27) to (33). Also, 

another zero located at s = -oo^ could be added to improve global 

stability as well as the speed of response of the simulator. 

If a zero is added to the modified analog equivalent position con- 
troller , i . e . , 


G f (s) - K r 
c 


(s/(Jj Q + 1) (s/to 1 + 1) (s/to 2 + 1) 
sCs/a)^ + 1) (s/w, + 1) 


(36) 


then the discrete-time equivalent position controller will still be 
given by equation (11), however, Kq, c Q , Cj, c 2> and c 3 will change to 
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K' T u) 0 m, 
3 4 


0 2 a)QW^a) 2 
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C 2 = 


>0-1) [(f) 2 - 
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(41) 


+ j (a) 3 +0)^) + 0)^0)^ 


The new values of the coefficients Kj to can again be computed 

off-line using equations (37) to (41) and (27) to (30). Notice that the 
values of the coefficients^ to Ky will not change because they do not 

depend on the zeros. It should also be pointed out that when a zero is 
added to the modified analog equivalent position controller the value of 
the coefficients through may be larger than 1 and can cause some 

overflow problems [5] if provisions are not made to handle such numbers 
because the INTEL 8086 microprocessor utilizes 2*s complement arithmetic. 
This problem can be alleviated by reducing the dc gain K' of the equiva- 
lent analog position controller of equation (36). 

Another possible improvement in the performance of the simulator 
could be achieved by introducing frequency prewarping of the analog 
equivalent position controller corner frequencies; however, the improve- 
ment would be minute in this case because the sampling rate used in the 
control algorithm is unusually high (100 Hz). 


Once the coefficients K 1 through K ? have been recomputed off-line, 

the old values can be changed interactively in the control computer 
when INIT is selected from the main menu. 

It should be apparent by now that the suggested changes only 
involve the recomputation of the coefficients Kj through K ? since they 

reflect any changes that may take place in the location of the poles and 
zeros of the modified analog equivalent position controller of equation 
(10) or the addition of another zero as in equation (36). 

Although the performance requirements of the simulator are not 
terribly stringent, it is deemed by the author that a much more robust 
controller can be designed to operate at a substantially lower sampling 
rate (much lower than 100 Hz as it is presently implemented) and yet 
outperform the present digital position controller in every aspect, such 
as global stability, speed of response, and flexibility. This, however, 
would require rewriting the control software. 
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To begin with, an improved controller design would have to take 
into account the dynamics of the resolver to digital converter which is 
used to determine the angular position of each of the simulator axes, 
since this device is installed in the feedback path of the position 
loop. This type of converter is modeled as a type II tracking loop and 
possesses the closed-loop transfer function (ILC Data Devices Corp.) 


g r/d 


(s) 


A(s/B + 1) 

1 3 2 A . . 

OB 3 + s + B S + A 


(42) 


where 

A = 150176 , XDC-19147-303 

| 184800 , XDC-1 9197-304 

B = \ 100 , XDC-19147-303 
| 300 , XDC-19197-304 


where the XDC-19147-303 and the XDC-19197-304 are the 14- and 16-bit 
resolution resolver to digital converters installed in the attitude 
simulator digital control cards, respectively. Furthermore, it would 
make use of the velocity signal which is also generated by the same 
resolver to digital converter since the noise content of this signal is 
much smaller than the one generated by a mechanical transducer such as 
a tachometer. 

This new controller would also implement the rate compensation 
totally in the digital domain since the velocity signal generated by 
the resolver to digital converter can be digitized using an analog to 
digital converter without much difficulty. This approach would add a 
tremendous flexibility to the controller design. 

Finally, this controller would be designed in such a way that its 
performance does not deteriorate with changes of simulator parameters 
such as inertia, i.e., it would have some level of adaptation that takes 
care of variations in the system parameters. 
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